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Abstract 

In this paper I develop a new computational method for pricing path dependent options. 
Using the path integral representation of the option price, I show that in general it is possible to 
perform analytically a partial averaging over the underlying risk-neutral diffusion process. This 
result greatly eases the computational burden placed on the subsequent numerical evaluation. 
For short-medium term options it leads to a general approximation formula that only requires 
the evaluation of a one dimensional integral. I illustrate the application of the method to Asian 
options and occupation time derivatives. 
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1 Introduction 



Financial derivatives (eg. options and futures) derive their value from an underlying traded finan- 
cial security, whose price is modeled by some stochastic process. In their most general form, the 
option payoff is path dependent since it depends on the entire future path traversed by the underly- 
ing security. Path dependent options are defined using either discrete or continuous price sampling. 
For continuous sampling closed from solutions are often available, but in practice most traded path 
dependent options are discretely sampled. It is known that the application of these closed form so- 
lutions leads to substantial pricing errors for discretely sampled options [jl], |2|, |3|]. This feature has 
necessitated the development of practical and efficient computational methods for the evaluation 
of path dependent options [Q]. Most research has focussed on either partial differential equation, 
Monte Carlo or tree based methods. In contrast to these approaches, in this paper I will develop an 
alternative based on the path integral formulation of the pricing problem. 

For many years, theoretical physicists have been developing and applying the path integral 
method for calculating expectations similar to those now being encountered in the evaluation of 
financial derivatives (via the risk- neutral valuation formula). A path integral is an infinite dimen- 
sional Riemannian integral as the integral is performed over a set of functions or paths. The path 
integral method is unique in that it gives a global formulation of the problem in question. This 
global description provides a powerful tool for deriving analytical approximation and numerical 
solution schemes that are difficult or impossible to formulate in other ways. Formally, the path 
integral method is easily extended to multi-dimensional problems. Much of the driving force be- 
hind the development of path integral numerical methods in physics has been that, compared to 
other methods, they show only a slow increase in computational complexity as the dimensionality 
of the problem increases. Path integral methods were first introduced by Feynman in 1942 as an 
alternative formulation of quantum physics They have found wide application in the evalu- 

ation of both the real time dynamics and equilibrium statistical mechanics of quantum many -body 
systems [0, @]. They are also a popular and natural tool for the analysis of diffusion processes 
[ |10| , |TTj , [I2p, including non-Markovian systems where there is a lack of practical alternatives [|14|]. 

The application of path integral methods to financial derivatives was pioneered by Dash who 
developed a path integral framework for pricing bonds and options within one factor term structure 
models [ |i5| , |T6| , [T7| ] . More recently, Linetsky JT8| ] was the first to show how path dependent options 
could be formally priced in a path integral framework. He considered several examples of one and 
two factor path dependent options. Baaquie JT^ ] has shown how path integral methods can be 
used to price vanilla options with stochastic volatility. The same author has also recast the popular 
Heath- Jarrow-Morton model of forward interest rates as a problem in path integration pOQ. Similar 
to Dash, Otto has more recently shown how to use path integration to price bonds and bond options 
under general short rate models [ pT] ] . Bennati et-al [ |2~2"| ] have focussed on a multi-dimensional path 
integral formalism for solving general financial problems based on systems of stochastic equations. 
All these preceding authors focussed on general formalism and exactly solvable models. They have 
shown that path integrals constitute a natural framework for describing the evaluation of general 
multi-factor derivatives. Although path integral methods offer an attractive way of obtaining exact 
solutions, they cannot find exact solutions not available using more standard methods. The greatest 
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promise of path integral methods will be in the development of new numerical and approximation 
methods for addressing pricing problems where exact solutions are impossible. 

Path integral numerical methods involve the evaluation of a multi-dimensional integral, either 
by deterministic or Monte Carlo methods. A review of deterministic and Monte Carlo methods, 
developed for physics applications, can be found in Drozdov [ |TT[ ] and Makri [[7]] respectively. An 
independent and much more recent development has been the application of path integral nu- 
merical methods to financial derivatives. This was pioneered by Eydeland [ J23] ] who, using fast 
Fourier transform methods, has derived a deterministic path integral algorithm for calculating the 
generating function of a random variable defined as the time integral of a general diffusion pro- 
cess. He pointed to a number of potential financial applications of this algorithm. Chiarella and 
El-Hassan have applied the method of Eydeland to the pricing of American bond options in the 
Heath-Jarrow-Morton framework p3]]. More recently, Chiarella et-al have devised a deterministic 
path integral algorithm based on Fourier-Hermite series expansions and applied it to the pricing 
of American options and point barrier options [ |Z5| , . They report computational times that are 
a significant improvement over the standard binomial and finite difference methods. In contrast 
to these deterministic methods, Makivic [ |2"%| ] has initiated research into the path integral Monte 
Carlo evaluation of financial derivatives. He broadly outlines a computational approach based on 
the standard Metropolis algorithm. Rosa-clot and Taddei [ |27| ] have discussed both a deterministic 
method and the Monte Carlo approach and applied them to several examples. Some of these pre- 
vious papers point out that the path integral method is superior to the the traditionally used lattice 
methods, because the underlying asset price is left continuous rather than being discretized. This 
has several important advantages. The option price is obtained more accurately since all possi- 
ble price paths are included in the simulation. Option Greeks are obtained more reliably since the 
method avoids the need for numerical differentiation. A further advantage of path integral methods 
is that they can be easily and efficiently extended to evaluate multi-factor financial derivatives. 

The application of path integral methods here is fundamentally different from all the previous 
cited works. We use the path integral representation of the option price to show that rather gen- 
erally, it is possible to perform analytically a partial averaging over the underlying risk-neutral 
diffusion process. This key result will greatly reduce the computational burden placed on the 
subsequent numerical evaluation. The application of this method is inspired by a technique first 
developed for computation problems in chemical physics [ J30] , [31] ]. Conceptually, the partial aver- 
aging corresponds to averaging over the high frequency fluctuations of the risk-neutral diffusion 
process. For short-medium term options it leads to a general approximation formula that only 
requires the evaluation of a one dimensional numerical integral. Longer term options can be evalu- 
ated deterministically or most generally by standard Monte Carlo methods. In this case, the partial 
averaging method will greatly reduce the required dimension of the Monte Carlo simulation. 

The outline of this paper is as follows. In section 2 we develop a formal path integral represen- 
tation for the evaluation of rather general path dependent options. In section 3 we show how the 
previous path integral can be numerically evaluated after performing analytically a partial averag- 
ing over the underlying risk- neutral diffusion process. In section 4 some examples are presented. 
For clarity of presentation and completeness, the path integral methods used in this paper are de- 
veloped from first principles in the three appendices. 
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2 Path dependent option theory 



In this section we will present a formal path integral framework general enough to price a wide 
range of path dependent options. We will assume the standard geometric Brownian motion (GBM) 
model of the asset price process. More general diffusion processes present no special difficulties. 
In this case the risk-neutral price process is given by the Ito stochastic differential equation 

dS t =(r-q)S t dt + aS t dW h (2.1) 

where r is the interest rate and q is the continuous dividend yield. Using Ito's lemma we can show 

dx t = fjdt + adW t , (2.2) 

where ^ 

x t = In S t , fi = r-q - -a 2 . (2.3) 

Consider a path dependent option with price Cp at expiry u given by 

C F (S u ,1,u) = F{S U ,1) = F(e Xu ,T), (2.4) 

where F is the option payoff function which depends on some path dependent random variable 1. 
In this paper we will assume that 1 can be written as 

1=1 dsw(s)f(x s ,s), (2.5) 



which is a time integral over an arbitrary function of the risk-neutral diffusion process (2.2). For 
continuous sampling w(s) = 1, but for discrete sampling (which is more realistic in practice) 

i 

where W{ are the sampling weights and Sj are the sampling times. The above definition of 1 
was used before by Wilmott et-al whose focus was on partial differential equation methods. 
It was shown to be general enough to include Asian, barrier and lookback options which are 3 
qualitatively different path dependent options. We will present examples in section 4. 
In a risk-neutral framework, the option price at inception time t is given by 

C F (S t ,t) = e-* r E Xt [F(e x °,T)], (2.7) 

where T = u — t and the expectation is with respect to the transformed risk-neutral price process 
(2.2) conditioned on the initial value x t . The price at inception can then be expressed as 

/oo roo 
dx u / dl P(x u ,1\x t )F(e Xu ,1), (2.8) 

where P(x u ,1\x t ) is the joint probability density function (PDF) of x u and the path dependent 
random variable 1. In appendix A we show for a general diffusion process how the joint PDF can 



4 



be formally computed as a path integral. For the special case of GBM, we show in appendix A that 
the joint PDF is given by 

/oo 
dke- ikX K(x u ,x t ;T), (2.9) 
-oo 

where fi is the constant defined in (2.3) and K, which we refer to as the propagator is defined by 

K(x u ,x t ;T) = Vx s exp 

Jxt 

In (2.10), the integration measure Vx s denotes a path integral which is defined precisely in ap- 
pendix A. It describes an infinite dimensional integral over all paths connecting x u at the expiry 
time and x t at the initial time. We refer to the function V in (2.10) as the potential function. For 
the GBM model it has the form 

V(x s , s) = -lika 2 w(s) f(x s , s). (2.1 1) 

We see that in this case the potential is imaginary with a functional form determined by the path de- 
pendent random variable (2.5). We show in appendix A that for more general risk-neutral diffusion 
processes the potential function is complex. 

The propagator (2.10), up to a boundary term, is the characteristic function of the joint PDF. 
In physics, the path integral (2.10) is equivalent to the path integral representation for the equilib- 
rium quantum statistical density matrix of a particle in a complex potential V(x) [Q], In this case 
temperature replaces the role of time in (2.10). Under an imaginary time transformation, (2.10) 
becomes equivalent to that which gives the quantum mechanical propagator of a one dimensional 
quantum particle in a complex potential V(x). These identifications with standard problems in 
theoretical physics are of great value because we can then use the methods and results of theo- 
retical physics for performing these path integrals. Tables of known exact path integrals of the 
form (2.10), for various potential functions, are listed by Grosche [^. Such tables along with the 
formulation provided here provide an easy way to obtain exact solutions to path dependent option 
prices. 

2.1 Seasoned path dependent options 

In this paper we consider the option at its inception time. More generally, the option price at time 
t' (£<£'< u) for a seasoned path dependent option is given by 

C F {S t ,,l(,t') = e- r ^E %t \F(e^X + J")], (2.12) 

where we use the more detailed notation 

Z u t = £ dsw(s)f(x a ,s). (2.13) 



P(x u , l\x t ) = 

Z7T 



l~tx u 



a*- 



y 2 T 
2o 2 



2o 2 



I ds (x 2 . + V(x s , s)) 



(2.10) 
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We then find 

, , , roc roo 

C F {S tl ,Tl,t') = ) / dx u / dZ$P(x u ,I$MF{e x ",I t t +1?,). (2.14) 

J — oo J — CO 

We see that this case only affects the payoff function in a simple way and the joint PDF we need 
to find is the same problem as before. Therefore all the final results can be trivially extended to 
seasoned options. 



3 Partial Averaging 

In the previous section we showed that for the GBM model, equations (2.5) and (2.9-11) define a 
formal path integral representation for the joint PDF. The option price is then obtained from this 
joint PDF via (2.8). In this section we use the previous path integral formulation to show that it is 
possible to perform analytically a partial averaging p3Dp over the underlying risk-neutral diffusion 
process. The option price can then be more efficiently evaluated by numerical methods. 

First we must discretize in time (2.2), by defining the discrete time s n = ne + t, where n = 
0, 1, N and e = T/N with T = u — t. The propagator (2.10) can then be decomposed as 



roo N 

K(x u ,x t ;T) = / dx N - 1 ...dxi'[[K(x n ,x n - 1 ;e) 

J -°° n=l 



(3.1) 



where x N = x u , x = x t and a general form for the short-time propagator is, as shown in appendix 
C, 



( 1 



V27ra% 



1/2 



cxp 



Xn-l) 



2(7% 



(3.2) 



where 7 will be defined below. Substituting (3.2) into (3.1) we find that the propagator becomes 



K[x u , x t ; T) 



1 



As e — > 0, its possible to show that 



dxN-i---dx\ exp 



JV 



-E 



n=l 



2(7% 



N 



l{x n ,x n ^i,e) 



n=l 



) + V(a; n _i, s„-i)) 
4a 2 



(3.3) 



(3.4) 



yields the correct short-time propagator when substituted into (3.2). 

We will refer to equation (3.2), with (3.4), as the primitive short-time propagator. It is the 
standard short-time propagator used in the numerical evaluation of path integrals. A key feature 
of this propagator is that it is not correct to first order in e. It is in fact only valid as e — > 0. 
Clearly, the dimension of the integral in (3.1) could be made much smaller by searching for short- 
time propagators accurate over larger time-steps. This observation has motivated the search for 
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improved short-time propagators for use in path integral calculations for physics applications [0, 
|13|, [3TJ]. In appendix C, we show that its possible in general to write 



00 1 / e \ m 

^fi-^ny X n — 1) ^) / . T ( ~Z n ) C m (x n , Ifi-li £J) 



(3.5) 



where C m describes a cumulant structure with each cumulant of order (o" 2 e) m_1 . The key result 
is, we can obtain a short-time propagator formally correct to second order in e, by truncating all 
cumulants beyond the first. This truncation corresponds to retaining an averaging over only the 
high frequency fluctuations of the risk- neutral diffusion process; i.e. a partial averaging. The 
improved short-time propagator will lead to a much more efficient numerical evaluation compared 
to that obtained by using the primitive short-time propagator. In appendix C, we calculate the first 
2 cumulants exactly for a general potential and derive an expansion of the propagator to third order 
in T (a result only valid for smooth potentials). 

Using the results from appendix C, with the GBM potential (2.11), we find that 



e 2 k 2 



where 



a(x r 



-ieka(x r , 



, Xr 



-i;e + 



■(5(x r 



drw( 



P(Pr) 



: exp 



,_i;e) + o(a e 
dp T P(p T )f(x T +p T ,r) 



and 



(3.6) 

(3.7) 
(3.8) 

t)t, x t = r(x n - x n -t) + x n -i, r = (s - s n -x)/£. (3.9) 

Equation (3.7) is a key equation as it contains the partial averaging. The origin of a is the first 
cumulant in the expansion (3.5) and its evaluation will lead to a short-time propagator correct to 
second order in e. Expanding a to order e is consistent with the truncation of the higher cumulants 
in (3.5). The origin of (3 is the second cumulant which is formally calculated in appendix C. All 
we need to know here is that (5 is of order a 2 e, since it will be set to zero at the end. We keep it to 
determine the order of the first correction term due to the truncation of all higher cumulants. After 
combining (3. 6), (3. 3) and (2.9) and performing the integration over k, we find that the joint PDF 
becomes 



a 2 e(l 



P(x u , 1\xt) 



dx 



N-l 



X 



exp 



dx x P[x Nj ..,xi\Xq} 

2^ 2 X/n=l P{%rn %n—l] ^) 



JV 



-1/2 



2lT£ 2 J2/3(Xm 



(3.10) 



where x N = x u , x t = x and the discrete path PDF is given by 



P[x N , ...,xib ] 



2ira 2 e 



exp 



N 

E 

n=l 



(Xr 



n=l 



2a 2 e 



(3.11) 
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Equation (3.11) describes the probability density of realizing a particular discrete path of the risk- 
neutral stochastic process (2.2). In the limit that j3 tends to zero, the joint PDF (3.10) becomes a 
delta function in T and we can show, using (2.8), that the option price becomes 

/■oo / N \ 

C F (S t ,t) ~ e~ rT / dx N ...dx 1 P[x N ,...,x 1 \x \ F[e XN , e ^ a{x n , x n -i, e) \ + o(s 2 a 2 T). 

J -°° n=l 

(3.12) 

The order of the correction term follows from (3.10) and a saddle point expansion of the Gaussian 
integral over J in (2.8). 

Equation (3.12) is the major result of this paper. The multi-dimensional integral can be evalu- 
ated by deterministic methods or more generally by standard Monte Carlo methods. In this case we 
use the observation that (3.12) is equivalent to an expectation of the option payoff function F, with 
respect to the discretely sampled risk-neutral diffusion process defined by (3.11), or equivalently 
by (2.2). In (3.12), the discretization time-scale e is independent of any discrete option sampling 
time-scale. If we choose e to match the interval between discrete option sampling, we find that 
(3.7) will reduce to a primitive short- time propagator and (3.12) will become 

,oo / N \ 

C F (S t ,t) = e~ rT \ dx N ...dx 1 P[x N ,...,xi\x ]F \e XN ) ^w n f(x n ) . (3.13) 

J ~°° \ n=0 / 

This is equivalent to a direct discretization of (2.7), consistent with a discretely sampled path 
dependent random variable described by (2.5) and (2.6). In this case no analytical partial averaging 
has been performed and the Monte Carlo evaluation of (3.13) is completely standard. The great 
advantage of the partial averaging method is that in (3.12), the discrete time interval e can be chosen 
to be much larger than the option sampling time-scale. The partial averaging is performed in (3.7) 
where we must average over Gaussian fluctuations about the straight line path x T connecting x n 
and x n -i. This corresponds to averaging over the high frequency fluctuations of the risk-neutral 
diffusion process. In practice, as will be seen in the next section, the partial averaging integral 
is easily performed analytically. It is simply the Gaussian transform of the function /, which 
defines the path dependent random variable in question via (2.5). Of most practical interest will 
be discretely sampled path dependent options, where w(s) is given by (2.6). In this case the 
subsequent integral over r in (3.7) reduces to a discrete sum which presents no problems. 

For the special case N — 1, for which e = T, x\ = x u and x = x t , we find that (3.12) 
becomes 

C F (S t ,t)^-^=J^dx u e W _ (^~^-^) 2 F ^^ Ta{XwXt . T ^ +o(a 2 T 3y 

(3.14) 

This describes rather generally an approximate path dependent option price. It can be simply 
evaluated as a one dimensional numerical integral. As a measure of the accuracy of (3.14), we ask 
at what time to maturity T, does the error term in (3.14) become 1% of the true option price. We 
assume a typical market volatility of a = 0.25 and that the unknown coefficient of the correction 
term is equal to the true option price. We find that T ~ 0.5 gives a 1% error, while T ~ 0.25 
gives an error of approximately 0.1%. These estimates begin to illustrate the power of the method 
presented here. 
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4 Examples 



The previous formulation is rather general and can be applied to a range of path dependent options. 
In this section we will show how two important and qualitatively different classes of path dependent 
options fit into this framework. 



4.1 Average rate options 

The payoff of the geometric Asian option is some function of the path dependent random variable 
given by 

ru 

1=1 dsw(s)x s , (4.1) 



where x s is related to the risk-neutral asset price by (2.3). From (2.5) we can identify / with x s . 
After performing the partial averaging (3.7), we find 

a(x n , x n -i, e) = [ drw(r)x T . (4.2) 
Jo 

For the continuous sampling (w(t) = 1) we find 

a(x n , x n _i; e) = (x n + z n _i)/2. (4.3) 

For this simple example a is just the primitive short-time propagator. 

The payoff of the arithmetic Asian option will be some function of the path dependent random 
variable 

ru 

dsw(s)e Xs . (4.4) 



From (2.5) we identify / with e Xs . After performing the partial averaging (3.7) we find 

a(x n ,x n - 1 ;e)= [ dr w(r)e^ +u ' /2 . (4.5) 
Jo 

We can expand (4.5) to order a 2 e without a significant loss of accuracy. This is consistent with the 
order of the truncation of the cumulant expansion (3.5). We then find 

a(x n , x n -i; e) ~ J dr w(r)e* T (l + a 2 e(l - r)r/2 + o(<7 4 £ 2 )) . (4.6) 

The final result will depend on whether we use discrete or continuous sampling. For continuous 
sampling we have w(r) = 1 and (4.6) becomes 

-(e a -l) +^(e>-2) + a + 2) + o(aV) , 



(4.7) 



where a = x n —x n -\. For discrete sampling we can perform the necessary summations analytically 
so the method is equally effective. It is instructive to compare (4.7) with the a that generates the 
primitive short-time propagator. This is obtained by approximating (3.7) by 

a(x n , x n _i; e) ~ ]- (f(x n ) + /(x n _i)) . (4.8) 



9 



One can see that (4.7) will only reduce to (4.8) in the limit a — > (when / = e Xs ). 

We can now use (3.12) to perform a Monte Carlo evaluation of the continuously sampled 
Asian option. Using (4.7) allows us to obtain accurate results by simulating only relatively low 
dimensional random paths of (2.2). This will avoid the problems associated with the Monte Carlo 
evaluation of these options [|53p. 



4.2 Occupation time derivatives 

A large and important class of path dependent options are those where the payoff depends on the 
time the asset price spends within a given region. This class of options have been referred to 
as occupation time derivatives and they have been discussed in detail by Linetsky p4j [JS] ] and 



Hugonnier Q3q]. The valuation of debt and contingent claims with default risk can also be placed 
in this class 

As a simple example, consider an option with a payoff that depends on the time spent below a 
possibly time dependent barrier B s . This time is a path dependent random variable which can be 
expressed as (and similarly for the above case) 

I = J U dsw(s)H(B s -e Xa ), (4.9) 

where H is a simple step function. We can therefore identify H with the function / in (2.5). After 
performing the partial averaging (3.7) we find 



f (*£»i) *^n— lj 



J dTw(T)N((lnB s -z T )/i/ T ), (4.10) 



where N is the cumulative normal distribution function. This function is difficult to integrate 
analytically. However the real case of practical interest is discrete sampling. In this case the 
integration reduces to a discrete summation and presents no problems. Barrier options are the 
most basic and well known example of options that fall under this framework. An interesting 
generalization is the step option which has a finite knock-out rate. This is motivated by risk- 
management arguments and a variety of possible payoff functions have been discussed [ fP] , fl5| ]. 
All variations are included in this framework since the payoff function is an arbitrary function of 
the path dependent random variable. 

Consider the example of derivatives depending on the occupation time between two barriers 
Bi and B 2 ^\ Included in this class are double barrier options p9] ] and the range products such as 



range notes and corridor options pQ , [41]]. This time is a path dependent random variable which 
can be expressed as 

J = J U ds w(s) [H(B 2 - e Xs ) - H{B X - e*')] . (4.1 1) 

From (2.5) we can identify the function / with the difference of two step functions. After perform- 
ing the partial averaging (3.7) we find 



•J 

J dr w{t)[N [{In B 2 -x T )/v T ) - N([\nB x - x T )/v T )\. (4.12) 



1 the occupation time outside these barriers can be trivially constructed from this time 
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Our framework is easily generalized to include derivatives dependent on several occupation times. 



5 Discussion and Conclusion 

The aim of this paper was to present a new approach to evaluating the price of path dependent 
options. We considered options with the general payoff function (2.4), contingent on a path de- 
pendent random variable expressible in the form (2.5). The key results of this paper were the 
general evaluation formula (3.12) and the short time to expiry approximation (3.14). They give the 
option price after performing analytically a partial averaging, defined by equations (3.7-9), over 
the underlying risk-neutral diffusion process (2.2). Since the method is analytically based, it also 
gives the order of the error made by choosing a finite time discretization. Specific examples were 
presented in section 4. 

In (3.12), the partial averaging allows one to choose the discrete time interval e to be much 
larger than the option sampling time-scale. We could, for example, imagine choosing e to be one 
month for an option with daily sampling. Clearly, the partial averaging method can greatly reduce 
the dimension of the integral in (3.12). This integral can be evaluated most generally by standard 
Monte Carlo methods. In this case the partial averaging will greatly reduce the dimension of the 
random paths to be simulated. In effect, it allows random simulations to be replaced with deter- 
ministic calculations. Standard methods to increase the efficiency of the Monte Carlo evaluation 
can still be used. These include variance reduction techniques the simulation of sample paths 
using the Brownian bridge process [ ft2"| ] or the use of quasi Monte Carlo sampling In- 
terestingly, quasi Monte Carlo sampling is known to be more advantageous for low dimensional 
numerical integrals. The partial averaging method can therefore increase the relative gains made 
by the implementation of quasi Monte Carlo methods. 



The framework presented here can be extended to multi-factor path dependent options Q18| , [22Q 
and to path dependent options dependent on several path dependent random variables. Recent 
work has shown that a path-independent option price, in a stochastic volatility environment, can 
be approximated by pricing a more complex path dependent option in the usual Black-Scholes 
framework. The volatility risk is included by an extra continuous path-dependent payout function 



which has the same form as (2.5) Q45| , |46| ]. This work suggests that the computational framework 
presented here might also be used to price path dependent options consistent with the market 
implied volatility. Further work is required in this direction. 

Path integral methods have long been developed and used as a computational tool in theoret- 
ical and chemical physics. Hopefully the work presented here will stimulate more interest in the 
application of these methods to problems in computational finance. 
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A Path Integral Representation 



Consider a stochastic process x s , which obeys the stochastic differential equation Q 

dx t = —g'(xt)dt + adW t . 
In our notation x s = x(s). We will rewrite this equation as 

it = -g'(xt) + o"Ct 



(A.l) 



(A.2) 



where Q = and a dot denotes the derivative with respect to time. This notation is more suited 
for the path integral formulation. In (A.l) we assume that o in independent of x t (additive noise). 
The path integral formulation of the more general case (multiplicative noise) can be found in fl44{]. 

Consider a general continuous Gaussian noise process x(s) (note that small s denotes the time 
history variable between time t and u and should not be confused with the price S t ). We discretize 
this process by defining a discrete time s n = ne + t, where n = 0,1, N and e = T/N with 
T = u — t. Note that e is equivalent to ds when N — > oo. The now discrete Gaussian process is 
fully defined by the normalized probability density functional 



Xn\ 



(2ttP /2 (deti?)~ i/z exp 



1/2 _ 



N 



^ Xn R n m Xr. 



n,m=l 



(A3) 



where \n = x( s n), E[x n Xm] = Rnm and denotes the inverse matrix. Equation (A. 3) fully 
defines the probability density of the whole history of the discrete Gaussian process. This normal- 
ization condition implies 

roc 

V [Xi, Xn] dxi-dxN = 1. (A.4) 



Since dW t has a variance dt, we know ( t will have a variance dt^ 1 . We then find that for the 
discrete time white noise process ( n , we have E[( n ( m ] = s~ l 8 nm where 5 nm is the unit diagonal 
matrix. Using (A. 3) we then find 



v[&,...,Cn] 

In the continuous limit this becomes 



N/2 



271 



N/2 



exp 



N 



Z n=l 



exp 



1 rU 



(A.5) 



(A.6) 



where the continuous time expressions are written with the understanding that e — > and N — > oo 
such that Ne = T. 

We wish to use the probability density functional (A.6) to find a probability density functional 
for the stochastic process x s . To do this we need to first discretize in time (A.2). In discrete time 
(A.2) becomes 

-g (x n ) + a( n , (A.7) 



•En— %n— 1 



e 



Note that any one-dimensional risk-neutral diffusion process can be cast into this form by a change of variable 
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where 

x n = 4>x n + (1 - 0)x„-i (A.8) 

and is a discretization parameter between and 1. From (A.7) we see that a path {d, ( N } 

maps to a unique path {xi, x^} as long as x is given. We will therefore write, by virtue of 
(A.4) 

/oo 

V[xi, xat|x ] J dxi...dx N = 1, (A.9) 



where J, defined by 



J = det 



9(n 



OXr 



m, n 



1,...,N 



(A.10) 



is the Jacobian of the change in coordinates. In the continuous limit we can substitute (A. 2) into 
(A. 6) to obtain 



V[x s 



■In 



N/2 



exp 



- 2 [ds (x s + 9 'i 



We show in appendix B that the Jacobian (A. 10) becomes in the continuous limit 



J 



I \-N 

[as) exp 



(f) ds g" (x s ) 



(A.ll) 



(A.12) 



The only sensible choice is = 1/2 [0]. Combining the Jacobian and (A.ll) we obtain a new 
probability density functional 



exp 



2a 2 



ds (x s + g'{x s ) 



2 1 

+ 2 



ds g"(x s 



(A. 13) 



which is normalized with respect to the functional measure V s which is the continuous limit of 



/ „\ -N/2 

T>x s = ylitso J dx\...dxN-\ 



(A. 14) 



This means the conditional probability density function (PDF) of (A. 2) is given by the path integral 



Dx s "P x \x s 



(A.15) 



The expectation of a functional J 7 fx,,], conditional on the initial value x t , can now be expressed in 
the path integral form 



E 



Xt 



F[x s 



dx u I Dx s x \P^s\ *^\p^s 



j-t 



Using (2.7), we can now write the option price in the path integral form 



C F (S t ,t) 



-rT 



dx u Vx s V x [x s }F(e x \l). 



(A.16) 



(A. 17) 



We will not deal directly with the above path integral representation of the option price. From (2.8) 
we see that we can extract the payoff function out of the path integral if we instead focus on the 
path integral representation of the joint PDF P(x u , T\x t ). In the next subsection we will see how to 
calculate this function. 
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A.l Calculating the joint PDF 

The joint PDF introduced in (2.8) can be simply expressed as the path integral 

P(x u ,l\x t ) = / Vx s V x [x s ]8(l - J). 



(A.18) 



Using the Fourier representation of the delta function and (2.5), we find that (A.18) becomes 

| roc 



P(x u ,l\x t ) = — dke~ lkI Vx s V x [x s ]exp 

Z7T J — oo J xt 

Substituting (A. 13) into (A. 19) and using 

ru 

J dsx s g'(x s ) = g(x u ) - g(x t ), 

we find that (A. 19) becomes 



dsw(s)f(x s , s) 



(A.19) 



(A.20) 



P(x u ,l\x t ) = — exp 

Z7T 



g(x t ) - g(x u ) 



/oo 
dke- ikI K(x u ,x t ;T), (A.21) 
-oo 



where K, which we will refer to as the propagator, is defined by 

1 



K(x, 



2a 2 



x s + V(x s , s) 



,x t ;T) = / Vx s exp 

Jxt 

and what we will call the potential function V is 

V(x s , s) = g' 2 (x s ) - a 2 g"(x s ) - 2ika 2 w(s) f(x s , s). 



(A.22) 



(A.23) 



Of special interest is the geometric Brownian motion model defined by (2.2). Comparing (2.2) 
with (A.l), we can identify g'(x t ) with /j, and we find the joint PDF (A.21) becomes (2.9) with the 
potential (2.1 1). Because the drift term is constant in this case, it can be extracted out of the path 
integral and we are left with an imaginary potential. 

B Calculating the Jacobian 

In this appendix we will calculate the Jacobian (A. 10). In discrete time the stochastic differential 
equation (SDE) (A. 2) becomes 



X n —X n —\ 



-g'{x n ) + a(x n )( n , 



(B.l) 



where it has been generalized to include a non-constant diffusion coefficient. From this discrete 
SDE we see that 

dC 

= 0, for m>n. (B.2) 
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Therefore the Jacobian matrix is triangular and from (A. 10) we obtain 



j - n dc " 



From (B.l) we find 



1 dx n 



where we have used 



0x n dx n dx n 



dX T 



obtained from (A.8). Multiplying (B.4) by ep 2 - we obtain 



For small e (B.6) becomes 

OXr 



d( n 8x n d( n d( n 
d( n _ 1 + £ <fig"(Xn) -£(j> (nCr'(Xn) 

dx n ea(x n ) 

~ ea(x n )(l -e<pg"(x n ) + e (j) Cnc'C^n)) 
~ £a(x n ) exp e <f>(a'(x n )( n - g"{x n )) 



Using this result we find 

N dx, 



N 



n 



n at 



n=l ^Cn \ n= l / 

In the continuous limit e — > 0, we therefore find 



AT 



n=l 



iV 



(j) j ds (cr'(x s )( s - g"(x s )^j 



J-'=e N i [ Y[a(x n )jexp 



From the continuous SDE (A. 2) we have 

Cs = (i s + g'{x s ))/a{x s ). 

Substituting this into (B.9) we find 



N 



-1 



J = e 



-JV 



n a &. 

\n=l 



exp 



(j) j ds {{x s + g'{x s ))a'{x s )/a{x s ) - g"{x s ) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



(B.ll) 



This agrees with equation 2.4.37 in Stratonovich [|10J who derives the Jacobian for a general multi- 
factor system of SDE's. For a constant diffusion coefficient (B.ll) reduces to (A. 12). 
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C Improved Short-time Propagator 



The goal of this appendix is to obtain an expansion of K(x u , x t ; T) in powers of T. The derivation 
here is inspired by the cumulant method used in connection with the quantum statistical density 
matrix [|3T]]. 

We must first write the propagator in the Fourier path integral representation (fpir). We de- 
compose the paths as 

/ 4a 2 T \ 1/2 oo 



where 



x T + 



7T 



E 

n=l 



Z„ Sill TLTCT 



n 



r(x u -x t )+x t , r = (s-t)/T. 



(C.l) 



(C2) 



The term x T in (C.l) is the straight line path connecting x t and x u . The remaining terms are the 
harmonic perturbations about the straight line path. With this we can now write 



2a 2 



ds 



x 2 s + V(x s ,s) 



(x u Xf) 

2a 2 T 



+ * E 4 + 



T r 1 



n=l 



2a 2 



dr V(x T , r). 



(C3) 



Summing over all paths between x t and x u is equivalent to integrating over all possible values of 
the Fourier coefficients {z n }. This means we can write 



/ — 

T>x s = (constant) x / dz\...dz 



(C4) 



where the constant is some Jacobian factor resulting from the change of integration variables. 
Substituting (C.4) and (C.3) into (2.10), we obtain the fpir of the propagator 



/oo 
dzi...dzooexp 
-oo 

where Kf(x u , x t \ T) is given by 

K f (x u ,x t ;T) -- 



-E 



T r 1 



n=l 



2a 2 Jo 



/ drV(x T . 
Jo 



T 



(C5) 



1 



2na 2 T 



1/2 



exp 



{x u — x t ) x 
2a 2 T 



(C6) 



The constant Jacobian factor in (C.4) must equal the prefactor of (C.6) to obtain the correct prop- 
agator when V = (clearly the Jacobian is independent of the potential). We can rewrite (C.5) 
as 



K(x u , x t ; T) = K f (x u , x t ] T) E 



exp 



T r 1 



2a 2 



/ drV{x T1 T 
Jo 



{Zn} 



(C7) 



where the expectation is with respect to the set of independent Gaussian random variables {z n } 
with zero mean and variance 1/27T. The FPIR can also be set up using a reference potential that is 
quadratic 
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We can use the FPIR (C.7) to derive an expansion in time for the short time propagator. We 
have using the cumulant expansion 

K(x u , x t ; T) = K f (x u , x t ; T) exp — (- —j C m (x u , x t ; T)j , (C.8) 



where 
and 



Ci = Mi, C, M 2 - M 2 . C 3 = M 3 - 3M 1 M 2 + 2Mf, 



Mm = E 



I cItV(x t ,t) 
Jo 



(C9) 



(CIO) 



{z n } 



The propagator (C.8) has the general structure defined by (3.2) and (3.5). Consider the first cumu- 
lant which is 



C x {x n , x t ; T) = £drE [V(x T , r)] {Zn} . 



(C.ll) 



We know from (C.l) that x T is a sum of Gaussian random variables which itself must be Gaussian 
distributed. So when calculating the expectation (C.ll), we can replace x T by 



where p T is a single Gaussian random variable with variance 

^2^g sm W)^ 2r(i _ T)T 



7T Z 



n=l 



(C.12) 



(C.13) 



Using this result we find that 

/oo 
^dp T P(p T )V(x T +p T ,r) 

where 



P(Pr) 



exp(-p 2 T /2u 2 T ) 



(C.14) 



(C.15) 



The first cumulant then becomes 

Ci(x u ,x t ;T) = dr dp T P(p T )V(x T + p T ,r). (C.16) 

JO J-00 

Expanding the potential around p T = and performing the Gaussian integrals we obtain 

2 4 

E [V(?r, r)] {Zn} ~ V(x T , t) + ^-V"(x T , t) + ^\/' w (a: r , r) + o(T 3 ), (C.17) 
which from (C.13) is an expansion in time T. 
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Consider now the second cumulant 



C 2 (x u , x t ; T) = jf 1 drdr' [e [V(x t , t)V(x t ,, r')] {Zn} - E [V(x T , r)] {Zn} E [V(x T ,, r% Zn} 



In (C.18) we can replace x T and x T > with 

x T = x T + p T , X T i = X T i + p T i 

where p T and p T > are two correlated Gaussian random variables with variances 

v 2 T = a 2 T(l - t)t, u 2 , = a 2 T(l-r')r'. 

Using (C.l) we find the co variance is 

,. ^ r , Aa 2 T ~ ^ r , sin(n7rr) sinfn'Trr') 

c(t,t') = E{p T p T ,] = —— E[z n z n ,] {Zn} ^ 



71 



n,n'=l 



nn' 



Using 

and the identity 
we find 



E[z n Z n >]{ Zn } = l^nn' 

2 ^2, sin(n7rr) sin(n7rr') 

2- ~2 = W ~ T 9h 



71 



n=l 



c(r,r')=a 2 Tr l (l-r g 



(C.18) 
(C19) 

(C.20) 

(C.21) 

(C.22) 
(C.23) 
(C.24) 



where 77 is the lesser of r and r' and r g is the greater of r and t'. Using (A. 3), we can write the 
probability distribution of p T and p T > as 



P(Pt,Pt>) = 



1 1 



: CXP 



v 2 ,p T - 1cp T p T i + z/ 2 p 2 , 
2(z/> 2 - c 2 ) 



(C.25) 



27r {Jy T , - c 2 
We therefore find that 

/oo 
dp T dp T/ P(p T , p T >)V (x T + p T ,T)V(x T > + p T > , t') . (C.26) 
-00 

Substituting this result and (C.14) into (C.18), we find that the second cumulant is 

C2(x u: x t ;T) = / drdr' / dp T dp T < P(p T ,p T >)V(x T + p T ,r)V(x T < + p T >,r') 

Jo \J-00 

dp T dp T , P( Pt )P( Pt/ )V(x t + p T , T)V(X T , + Pr ,, t')) . (C.27) 

Expanding the second cumulant around p = 0we find that the first order remaining term is 

C 2 (x u ,x t ;T) ~ a 2 T C drdr' V'(x T , t)V'(x t ,, t') r,(l - r g ) + o(a 4 T 2 ). (C.28) 
Jo 
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This follows since the integral we have to do is just the covariance c. 

It can be shown that the mth cumulant is of order (a 2 T) m ~ l . Therefore from (C.8) we see that 
C2 starts contributing at order T 3 , while C3 starts contributing at order T 5 . This means that to get 
the expansion (C.8) correct to order T 3 , we need only expand C\ to order T 2 and C2 to order T as 
has been done in (C.17) and (C.28). This means that C\ alone will give the correct propagator to 
order T 2 . Substituting (C.28) and (C.17) into (C.8) we find 



K(x u ,x t ;T) ~ 



K f (x u ,x t ;T)exp 
a 2 T 5 



T 



2a 2 



f 1 dr V(x T , r) - — I 1 dr r(l - r) V"(x T , r) 
jo 4 Jo 



+ 



16 

8^ 



f 1 drr 2 (l-r) 2 V" f '(x T ,r) 
Jo 

f 1 drdr' V\x T , t)V'(x t/ ,t') r,(l - r fl ) + o(T 4 ) 

JO 



(C.29) 



which is our desired expansion of the propagator to third order in time. 
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